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We show that the Almeida-Thouless line in Ising spin glasses vanishes when their dimension 
d -¥ 6+ as h AT /T* = C(d - 6) 4 (1 - T/r c ) d/2_1 , where C is a constant of order unity. It is shown 
that replica symmetry breaking also stops as d — > 6 + . Equivalent results that could be checked by 
simulations are given for the one-dimensional Ising spin glass with long-range interactions. 
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I. INTRODUCTION 



close to T c (0), as 



The field of spin glasses is now well into its fourth 
decade, with many important questions still unresolved. 
Chief among these is the nature of the ordered phase. 
The principal rival theories of it are (i) the replica- 
symmetry-breaking (RSB) theory of Parish, motivated 
by the exact solution of the Sherrington-Kirkpatrick (SK) 
mean- field model, and (ii) the droplet /scaling theorjPEI 
based on the properties of excitations in the ordered 
phase. There is still no general consensus on which (if 
any) of these theories is correct. 

An important (and perhaps the simplest) discrimina- 
tor between the theories is the predicted behavior of 
the system when the temperature is decreased in the 
presence of an applied magnetic field. According to 
the RSB scenario, there will still be a phase transition 
in an applied field h, ocurring at a temperature T C (K), 
the Almeida-Thouless (AT) line^, which decreases as the 
field increases. This line can be explicitly calculated in 
the mean-field SK model. For small h it has the form 
T c (h) = T c (0) - const. h 2 / 3 or, equivalcntly, 



If 



(1) 



(T c (0) = T c ). In the RSB theory, this line is the phase 
boundary between the (high-temperature) replica sym- 
metric phase and the (low-temperature) broken replica 
symmetry spin glass phase. Within the droplet/scaling 
theory, on the other hand, there is no such line: An 
applied magnetic field is predicted to remove the phase 
transition completely and the low-temperature phase in 
the absence of a field has replica symmetry. 

Before presenting our calculation, we recall that the 
shape of the AT line starts to differ from Eq. ( I]) already 
for d < 8, as shown by Green et alP and by Fisher and 
SompolinskjP. Instead of Eq. ([I]), these authors show 
that the AT line has the form 



h AT 



_ rpld/2-l 



[T c - T] 



for 6 < d < 8. Note that as d -t 8" in Eq. Q, Eq. Q is 
recovered. 

In this paper we will derive an exact result for the form 
of the AT line in d — 6 + e dimensions, correct to leading 
non-trivial order in e. The result can be written, for T 
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= C(d-6f 1 



T 



d/2-l 
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where C is a non-universal constant. The form is the 
same as that proposed by Green at al. and by Fisher 
and Sompolinsky, but the amplitude contains the factor 
(d — 6) 4 which vanishes (rapidly) for d — >■ 6. 

In itself this result does not prove that there is no 
AT line below six dimensions. If there were an AT line 
below six dimensions, it would have the scaling form 
h 2 AT ~ (1 — T/T c )^ +1 , where (3 and 7 are the critical 
exponents of the zero-field spin glasiP. The vanishing of 
the amplitude in Eq. |3]) would arise to ensure continu- 
ity in six dimensions of the forms above and below T<Pl 
Hence we need to provide an additional argument why 
there should be no AT line when d < 6. 

Even before the droplet /scaling theory had been devel- 
oped, Bray and Roberts^ (BR) had used standard renor- 
malization group (RG) methods to investigate the pu- 
tative RG fixed point that controls the critical behavior 
at the AT line. In zero field there are three degener- 
ate soft modes at the critical point (usually called the 
"longitudinal", "anomalous" and "replicon" modes). A 
conventional RG analysis- shows that the upper critical 
dimension is d u — 6. For d > 6 the Gaussian fixed point 
is stable and the critical exponents take their mean-field 
values. For d < 6, a non-trivial fixed point is stable and 
non-classical exponents are obtained and can be calcu- 
lated as a power series in (6 — d) in the conventional way. 

An applied magnetic field, however, changes every- 
thing. The field lifts the degeneracy, leaving a single soft 
mode, the replicon mode. In their RG calculation, BR 
discarded the two hard modes, and derived RG recursion 
relations appropriate to the soft modes which we discuss 
below. Again, the critical dimension is d u = 6. There are 
now two coupling constants Wi and ui 2 . In dimensions 
d > 6, the RG flows have a single stable fixed point, the 
Gaussian fixed point with w* = = 0, implying that 
for d > 6 there is a continuous phase transition, with 
critical exponents given by their mean-field values. 

For d < 6, however, no physical stable points could be 
found. BR suggested that this might be because there 
was no AT line for d < 6. We shall strengthen this ar- 
gument by examining the RG flows when d > 6 in more 
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detail. We find that the basin of attraction of the Gaus- 
sian fixed point shrinks to zero as e —¥ 0; it has a linear 
extent of order e 1 / 2 . It thus seems natural to expect that 
there is no AT line when d < 6: There is no physical 
stable fixed point when d < 6 and the basin of attrac- 
tion for the controlling fixed point in 6 + e dimensions is 
shrinking to zero as e — > 0. 

Above the AT line the high-temperature phase is 
replica symmetric. Below it, the phase has broken replica 
symmetry. Thus if the AT line is disappearing as d — > 6 + , 
then one would naturally expect that replica symmetry 
breaking in the zero-field case would vanish as d — > 6 + . 
We shall demonstrate that is indeed the case by show- 
ing that the "breakpoint" x\ in Parisi's RSB function 
g(x^goes to zero in this limit. Thus the low-temperature 
phase for d < 6 should be replica symmetric. 

The structure of the paper is as follows. In section II, 
we present an RG analysis, valid for d > 6, that leads 
directly to Eq. In section III we analyse the conse- 
quences of our RG results for the form of the breakpoint, 
x\ , in the Parisi RSB theory. We find the x\ tends to zero 
for d — > 6 + , suggesting that replica symmetry breaking 
goes away in six dimensions. In section IV, the RG equa- 
tions of BR are presented. The basin of attraction of the 
Gaussian fixed point is determined numerically. It is a 
compact region enclosing the origin, with linear dimen- 
sions of order y/e, shrinking to a point as d — > 6 + . In 
section V, we treat the one-dimensional spin-glass with 
long-range inetractions falling off with distance r as r~ a . 
This model has mean-field behavior for a < 2/3 and non- 
classical behavior for a > 2/3. Thus a < 2/3 corresponds 
to d > 6 in the short-ranged model. We show that the 
de Almeida-Thouless line goes away as a — ¥ 2/3~, and 
that the breakpoint X\ of the Parisi function goes to zero 
in this limit. We conclude with a brief discussion of the 
main points in section VI. 



II. RENORMALIZATION GROUP ANALYSIS 

We now describe the calculations that lead to the re- 
sult quoted above.. We start from the Ginzburg-Landau- 
Wilson free-energy functional for the Ising spin glass 
which, written in terms of the replica order parameter 
field, is 



F[{Q a p}\ = J d d x [\r J2 a<0 Q 2 ap + \ £ Q< ^(VQ Q/3 ) 2 

+ 0(Q 4 ) (4) 



where h is the applied field. We use conventional RG 
methodd^, but work above the upper critical dimension, 
d u = 6, so we define e = d — 6. A simple scaling analysis 
of the terms in the functional of Eq. Q shows that the 



use of the simple RG equations of the h = theory to 
obtain the AT line as T -> T c . When d < 6 this will 
not be possible and the full set of RG equations in Ref P3 
would have to be solved instead. 

The renormalization group flows for w and r read^, 



dw 1 r . ,, o 



rj(l)]r — Aw 2 r, 
while h 2 evolves according to 



d{h 2 ) _ d+2-r](l) 
dl ~ 2 



(5) 
(6) 

(7) 



where rj(l) = ~(2/3)w(l) 2 . The RG equations are correct 
to 0(w 3 ), and to linear order in r. There is an additional 
r-independent term of order w 2 in Eq. |6| that we have 
omitted since it ultimately just leads to a shift in r, i.e. a 
shift in the critical temperature. In these equations the 
usual geometrical factor = 2/r(d/2)(47r) d / 2 has been 
absorbed into w 2 . For d > 6 the Gaussian fixed point 
is stable and the critical exponents take their mean-field 
values. 

Integrating the RG flow equations when d > 6 up to 
scale I gives 



Mi) 

r(l) 

Hi) 2 



w(0)exp(-el/2) 



[(2w(0) 2 /e)(l-e- d ) + l] 1/2 

r(0) exp[2Z - (10/3)A(/)] , 

/i(0) 2 exp[{(d + 2)/2}Z + (1/3)A(0] 



(8) 

(9) 
(10) 



where 



A(0 = / w{l'Ydr = - In 
Jo 2 



and to(0) = w, h 2 (Q) = h 2 , r(0) = r. 
At large I, 



w(l) 



2(1 + e/2w(0) 2 ) 



1/2 



cxp(-d/2). 



(11) 
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The exponential follows from the 0(w) terms in the RG 
equation jsl. The 0(w 3 ) terms in Eq. ^ serve to fix 
the amplitude of the asymptotic decay. The basic idea 
is to integrate the RG equations to a scale I* at which 
the running coupling constant, w(l*), is small enough for 
one-loop order perturbation theory to be accurate. 

The one-loop perturbative calculation of the AT line 
has been carried out by Green et al.. The result id^ 



h 2 /Q = lUw 



2, r\ 2 I d 



(13) 



natural size of h 2 is - \r\ 2 /w\ this remains the correct where I d is the integral J °° d d qq A {q 2 + |r|) 2 , which 
scaling form for all d > 6. h 2 AT is given by Eq. ^ and equals A. 



d\r 



d/2-4 



and A d = 1/2 for d = 6. In Eq. (|13j, 

Eq. ( 16 1 below, and is much smaller than |r| 2 /w provided the geometrical factor K d has again been absorbed into 

W' 



\ryT^ <C 1. We shall work in this limit as it allows the w 2 , and Q is the mean-field value of the order parameter 
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for r < 0. The factor of 36 between this result and that 
irP is due to the different definitions of w, which differ by 
a factor of 6 in the two calculations. Using the result Q ~ 
3|r|/iu, which holds on the AT line in the perturbative 



limit where 



K 2 « 1, we find 



w\r\ 



d/2-l 



(6<d< 



(14) 



which is the result obtained in RefsP ancP. 

To exploit the perturbation theory result, we coarse- 
grain to scale I and replace w by w(l), r by r(l) and h 2 by 
h{l) 2 to obtain h{l) 2 ~ w(l) |r(0| d/2_1 . Eq. O shows 



that w(l) becomes progressively smaller as I increases. 
Our use of the perturbative result will become valid at 
some value 1 = 1*, where I* will be specified below. In- 
serting the I dependent forms for h{l), w(l) and r(l) at 
1 = 1*, one obtains 



MO) 5 



w(0)|r(0)| 



d/2-1 



L AT 



[(2 W (0) 2 /e)(l-e- d *) + l] 



5(2/6-1 



(15) 



Setting d = 6 + e, and taking exp(— el*) <C 1 we obtain 
the final result, correct to leading order in e, that 



l AT 



w(0)\r(0)\ 



d/2-l _ ( 6 

" V2w 2 



This is equivalent to Eq. ([3]). 



\d/2-l 



(16) 



One cannot simply set I* — > oo in Eq. (15). The 



RG equations as presented here are valid provided r(l*) 
remains small or comparable to the square of the cut-off 
in the theory, which is conventionally taken to be unity. 
To leading order r(l*) = r(0) exp(2Z*), so we shall fix I* 
by setting^ 



\r(l*)\ « |r(0)|e 



21* 



1. 



Eq. ( 16 ) will hold provided that 



"(0)1 



:/2 = 



|r| £ / 2 « 1. 



(17) 



(18) 



Thus as e — > 0, the temperature interval near T c over 
which Eq. ([3| is accurate becomes very narrow. 

III. DOES RSB DISAPPEAR FOR d < 6? 

We are arguing in this paper that the lower critical 
dimension for the AT line - the line beneath which RSB 
sets in in the presence of a field - is 6. In zero field there 
is still a phase transition when d < 6, provided d > di, 
where di is the lower critical dimension, thought to be 
between 2 and 3. One might wonder whether replica 
symmetry breaking is present in the zero-field phase or 
whether RSB too disappears below six dimensions. We 
next present an argument that RSB in zero field vanishes 
as d — > 6 + based on the same methods which we used to 
derive the dimensionality dependence of the AT line. 

When d > 6 the form of replica symmetry breaking 
is well-establishecP^. {Q a p) becomes a function q(x) in 



the interval 1 > x > , which is constant for 1 > x > x\ 
at the value qea (= Q) and then falls from this value to 
zero at x = 0. For the functional of Eq. |4) (with h 2 set 
to zero) the breakpoint x\ is given to one- oop order for 
6 < d < 8 bjffH 

x x - w 2 |rf/ 2 - 3 . (19) 

If x\ were zero, the spin glass would be replica symmetric. 
We shall argue that as d — > 6 + , x\ goes to zero linearly 
with (d— 6), suggesting that when d < 6 there will be no 
replica symmetry breaking. Higher loop terms will leave 
the exponent of \r\ unchanged, but modify its prefactor, 
just as they do for the AT line of Eq. 

We shall work in 6 + e dimensions again with e small 
and take |r| small so as to permit neglect of the quartic 
terms in the functional Eq. Q. The RG will be used 
to coarse-grain to scale I* at which the form of the per- 
turbative one-loop order expression for the breakpoint X\ 
becomes valid. The breakpoint x\ [r, w] has zero scaling 
dimension, so it follows that 



Xl [r, W ] = xMnMn] ~ Hn 2 Hn\ d/2 ~ 



(20) 



Inserting our previous expressions for r(l*) and w(l* 
one finds 



Xi\r,w\ 



w{Q) 2 \r(0)\ d / 2 - 3 
[(2w(0) 2 /e)(l _ e -d-) + i]5d/6-4- 



(21) 



I* is specified as in Eq. (18). Hence in the limit when e 
is small (but el* > 1) 



Xi 



hd-&)\r\ d ' 2 -\ 



(22) 



which goes to zero as d — > 6 + , implying that replica 
symmetry breaking vanishes in this limit . 



IV. THE BRAY-ROBERTS CALCULATION 
AND THE ISLAND OF STABILITY 

We shall now describe the BR calculation of the RG 
equations pertinent to the AT line. It is this work which 
is at the heart of our contention that 6 is the lower critical 
dimension for the existence of the AT line. As mentioned 
before, their approach is to study just the fields in the 
rcplicon sector Qa/3, which are such that ~^2pQ a p = 0. 
The effective functional is 



F[{Q aP }} = f d d x 



+ (fl/6) QapQp-fQ-fa + (w 2 /6) Qlf 



(23) 



Here the convention has been adopted that the sums over 
replica indices are unrestricted except that Q aa = 0. At 
the AT line, f = in the mean-field approximation. Ac- 
cording to BR, the coupling constants W\ and W2 would 
have for d > 6 the RG flow equations 

\ [-£ - 3t7k(0] wi + Uwl-36w 2 lW 2- 



(II 



-I8wiw 2 +W2 
(24) 
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FIG. 1: The island of stability of the Gaussian fixed point 
for d > 6. Distances are measured in units of \fk. Only the 
region W2 > is displayed. 



-e — 3rjn(l)] W2 + 24wfw2 — GOwiw^ + 34w 



d,W2 __ 

~dT ~~ _ 

(25) 

Once again we have adopted the convention of absorb- 
ing the geometric factor Kd into w\ and u>2 and t]r(1) = 
(Aw\ — 1&W\W2 + llwl)/3. Presumably one could ob- 
tain these replicon sector equations by integrating the 
full equations of RefPS containing the hard longitudinal 
and anomalous modes as well as the replicon modes until 
the hard modes are decoupled from those in the replicon 
sector. If one were able to carry out this formidable task 
the initial values of w\ and W2 could be determined. We 
suspect that the initial value of w\ would turn out to be of 



order w(l*)(^ 



cl/2| r |e/4^ 



At mean- field level W2 is ~ yQ, 



where y is the coefficient of a particular quartic term, 
J2a<p Qt/3' i0un d when one goes beyond the cubic 
functional in Eq. (|4pl The effective value of y is domi- 
nated at small \r\ when 6 < d < 8 by its renormalization 
by the four cubic vertex "box" diagrarrP. Then y ~ 
w 4\ r \d/2-4 and so W2 „ «;(f*) 3 |r(i*)] 6 / a (~ e 3 / 2 |r| 3e / 4 ). 
Fortunately we do not need the precise initial values of 
wi and W2 for our argument, which is based upon the 
form of the basin of attraction of the Gaussian fixed point 
ofEqs. p4|-p5l). 



The basin of attraction has been determined numeri- 
cally, and is displayed in Figure 1. By scaling both w\ 
and W2 by e 1 / 2 and writing I — el, the explicit depen- 
dence on e in Eqs. ( 24 )- ( 25 1 can be removed. There is a 



peculiarly shaped compact region (the "island of stabil- 
ity" ) around the origin in the W\ , W2 space inside which 
all flows are to the Gaussian fixed point w* = = w^. 
Outside this region the flows are to infinity. The lin- 
ear extent of this basin of attraction is from the scal- 
ing of wi and W2 of order e 1 / 2 . Thus as e — > 0, the 
size of the basin of attraction shrinks to zero. Now BR 
showed that when d < 6 there was no stable physical 
fixed point; all flows of W\ and W2 were to infinity. Since 
the basin of attraction of the fixed point associated with 
the AT line is shrinking to zero as d —¥ 6 + and no phys- 
ical stable fixed points exists for d < 6, it seems natural 



to expect that the lower critical dimension of the AT 
line and replica symmetry breaking generally must be 6. 
This is a result in line with earlier expectations^ but is 
also in accord with more recent results for the "strongly- 
disordered spin-glass model", a model with an unusual 
bond-distributiorf^. 

The behavior of the basin of attraction of the Gaus- 
sian fixed point of Eqs. (24)-(25) as d — >• 6 + should be 
contrasted with that of Eq. (|5) for the zero-field case. In 
the latter case the basin does not shrink to zero as d —> 6 + 
and there is a phase transition when di < d < 6. 

Our treatment of the AT line was done using the zero- 
field RG equations with one coupling constant w rather 
than the full set of coupling constants, W\ , w% , • • • , w& of 
RefP^l. The justification for this is that all but wi can 
be dropped when obtaining the form of the AT line as 
T — > T c . To elaborate this point further we have already 
quoted forms for W2H*) and wi(l*) and as \r\ —> 0, 11)2(1*) 
is indeed ncglible in comparison with wi(l*). However 
for d < 6, the situation is quite different. If the coupling 
coefficients W2, W3, ■ ■ ■ ,Wg are anything but zero, the RG 
flows will take them to infinity and so for d < 6 it is never 
possible to work with just the one-coupling constant RG 
equations. The same point is relevant for the attempts 
to extend RSB calculations to d < 6 in zero-fielcPR For 
d > 6 all the coupling constants can flow to the Gaussian 
fixed point. 



V. THE ONE-DIMENSIONAL SPIN GLASS 
WITH LONG-RANGED INTERACTIONS 

Because our calculation on the form of the AT line 
is only valid when e is small, there is no chance that it 
can be directly checked by simulations. However, we can 
derive the analogous results for the one-dimensional Ising 
spin glass with long-range interactions, whose AT line 
has recentl y bee n the subject of contradictory numerical 
simulation d 13 ! 14 !. The Hamiltonian of these studies are 
variants of 



H 



<ij> 



ijSiSj 



(26) 



where hi is a random field of variance h 2 , the sum is 
over all pairs < ij >, and i and j are positions on the 
one-dimensional lattice. The interaction 



J i 1 J 1 



-J\ 



(27) 



where the ey are independent random variables with 
a Gaussian distribution of zero mean and unit vari- 
ance. This model was introduced by Kotliar et alP^, 
who showed that for a < 2/3 the model has mean- 
field critical exponents, and non-mean field exponents for 
2/3 < a < 1. When a > 1, there is no finite temperature 
phase transition. The a- interval, 2/3 < a < 1, is the ana- 
logue for short-range spin glasses of the dimension range 
between the upper critical dimension (d u = 6) and the 
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lower critical dimension, while a < 2/3 corresponds to 
dimensionalities d > 6. Thus just by changing the value 
of a one can explore both systems corresponding to high 
and low dimensionality. 

Our expectation is that when a < 2/3 there will be 
an AT line and both R e f s .lMiH confirm this. Unfortu- 
nately for a > 2/3 the two groups of simulators were 
in disagreement with each other: Refp21 did not see an 
AT transition, whereas RefP^I did. Our prediction below 
of the form of the AT line as a — > 2/3~ supports the 
conclusions of Refp21. 

The RG equations near the upper critical value of <r, 
2/3, were first written down in RefPP They exploited 
the fact that for these long-range interactions 2 — rj = 
2a- 1 16 . The bare propagator is l/(q 2<T_1 +r). The RG 
flow equations become 



dw 
~dl 
dr 
~dl 

while h 2 grows as 



-(2 - 3a)w - 2w 3 
(2a - l)r - 4w 2 r, 



d(h 2 
~dT 



ah 1 



(28) 
(29) 

(30) 



Eq. (30) is exact. 



Eqs. (28) and (29) are valid only when w is small, but 



The perturbative calculation along the lines of Green 
et alP of the AT line to one loop order gives a result 
similar to Eq. (14): 



h 2 /Q 



\^w 2 \r\ 2 l n 



(31) 



where I a = J °° dqq 2ly2<J 1 ^(q 2a 1 + \r\) 2 , which equals 
Bff | r |(5-8o-)/(2<7-i) ) and £ 2/3 = 3 _ Hence to one-loop 
order the equation of the AT line is 



i i 2 ~ 2 ° 
w\r \ 2a - 1 . 



(2/3 > a > 5/8). 



(32) 



For a — > 5/8, the mean-field AT form h 2 oc |r| 3 is recov- 
ered. Hence a = 5/8 is the analogue of 8 dimensions for 
short-range spin glasses. 

It is straightforward to integrate the RG equations for 
the long-range case to obtain the form of the AT line as 
a -> 2/3". The result is 



m ~ (w) 

or equivalently 

C(2-3cr) 5 ^ ?r ) 



iu(0)|r(0)| 



2 2:r 

2(7 — 1 



h 2 
T? 



T 

Tr. 



2 ■ 2 :T 

2a-l 



(33) 
(34) 

and C is again a constant of 0(1). Hence as a —> 2/3" 
the AT line goes away. That there was no AT transi- 
tion in the interval 1 > a > 2/3 has also recently been 



suggested^ from an expansion about a = 1, the "lower 
critical value" . 

The equivalent result for the "break-point", x\, in the 
Parisi function for the one-dimensional long-ranged sys- 
tem is 



xi ^(2-3a)|r|( 4 - 6ff )/( 2 - 1 ), 
when a — > 2/3—. 



(35) 



VI. CONCLUSION 

We have presented arguments that the Almeida- 
Thoulcss line in spin glasses is absent in systems with 
six or fewer space dimensions, i.e. these systems exhibit 
no phase transition under cooling if an external mag- 
netic field is present. We have also argued that the fea- 
tures associated with broken replica symmetry, in the 
Parisi solution of the SK model, are not present in finite- 
dimensional spin glasses with d < 6. 

Equivalent results have been obtained for a one- 
dimensional spin glass with interactions decaying with 
distance r as l/r CT . These systems are more amenable 
to simulation than short-ranged systems in high- 
dimensional space. To date, however, there seem to be 
no simulational studies of the one-dimensional long-range 
model which might help to confirm Eq. ( 34 ) in the inter- 



val 5/8 < c < 2/3. One issue which will complicate such 
studies was pointed out in Refp2J. Simulations at the AT 
line require the system to be large enough so that the 
Parisi overlap function P(q) has only positive support. 
This requires h(QN) 1 / 2 > T. Whether simulations can 
be done with the number of spins N large enough to meet 
this requirement remains to be seen. 

The debate as to the nature of the spin glass phase in 
three dimensional systems has run for so long because 
on the experimental side, for ex amp le, dynamical effects 
can produce an apparent AT lin o 3 ! 18 !, while in simulations 
there are always finite size effects which can mimic some 
of the effects of replica symmetry breaking^!. As a con- 
sequence it was always hard to be certain which of the 
two pictures, RSB or droplet scaling, was correct. Ref.^ 
is an example of a recent simulation and provides fur- 
ther references. We believe that our calculations provide 
strong arguments that the spin glass phase will not have 
replica symmetry breaking in dimensions d < 6. 
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